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Abstract 

The temperature dependence of the chiral condensate in isospin-symmetric nuclear matter at 
varying baryon density is investigated using thermal in-medium chiral effective field theory. 
This framework provides a realistic approach to the thermodynamics of the correlated nuclear 
many-body system and permits calculating systematically the pion-mass dependence of the free 
energy per particle. One- and two-pion exchange processes, A(1232)-isobar excitations, Pauli 
blocking corrections and three-body correlations are treated up to and including three loops 
in the expansion of the free energy density. It is found that nuclear matter remains in the 
Nambu-Goldstone phase with spontaneously broken chiral symmetry in the temperature range 
T < 100 MeV and at baryon densities at least up to about twice the density of normal nuclear 
matter, 2p ~ 0.3 fm~ 3 . Effects of the nuclear liquid-gas phase transition on the chiral condensate 
at low temperatures are also discussed. 



The chiral condensate (qq), i.e. the expectation value of the scalar quark density, plays 
a fundamental role as an order parameter of spontaneously broken chiral symmetry 
in the hadronic low-energy phase of QCD. The variation of (qq) with temperature and 
baryon density is a key issue for locating the chiral transition boundary in the QCD phase 
diagram. The melting of the condensate at high temperatures and/or densities determines 
the crossover from the Nambu-Goldstone phase to the Wigner-Weyl realization of chiral 
symmetry in QCD. 

It is thus of principal interest to perform a systematically organized calculation of the 
thermodynamics of the chiral condensate. Such a calculation requires knowledge of the 
dependence of the free energy density on the light quark mass (or equivalently, on the 
pion mass). The appropriate framework for such a task is in-medium chiral effective held 
theory with its explicit access to one- and two-pion exchange dynamics and resulting 
two- and three-body correlations in the presence of a nuclear medium. 

Previous studies of the in-medium variation of the chiral condensate were mostly 
concerned with the density dependence of (qq) at zero temperature, using different 
approaches such as QCD sum rules [1] or models [2,3] based on the boson exchange 
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phenomenology of nuclear forces. Temperature effects have been included in schematic 
Nambu-Jona-Lasinio (NJL) approaches [4,5]. Such NJL models work with quarks as 
quasiparticles and provide useful insights into dynamical mechanisms behind sponta- 
neous chiral symmetry breaking and restoration, but they do not properly account for 
nucleons and their correlations, a prerequisite for a more realistic treatment. 

The present work extends a previous chiral effective field theory calculation [6] of the 
density-dependent in-medium condensate (qq)(p) to finite temperatures T. Corrections to 
the linear density approximation are obtained by differentiating the interaction parts of 
the free energy density of isospin-symmetric nuclear matter with respect to the (squared) 
pion mass. Effects from one-pion exchange (with m^-dependent vertex corrections), it- 
erated l7r-exchange, and irreducible 27r-exchange including intermediate A(l 232)-isobar 
excitations, with Pauli-blocking corrections are systematically treated up to three-loop 
order. The dominant nuclear matter effects on the dropping condensate are supplemented 
by a further small reduction due to interacting thermal pions. To anticipate the result: 
we find that the delayed tendency towards chiral symmetry restoration with increasing 
baryon density p, observed at T = [6] in the same framework, gets gradually softened 
with increasing temperature. An approximately linear decrease of the quark condensate 
with increasing p is recovered at temperatures around T ~ 100 MeV. However, no rapid 
drive towards a first order chiral phase transition is seen, at least up to p < 2 po where 
po = 0.16 fm -3 is the density of normal nuclear matter. The only phase transition known 
for nuclear matter in this density range is the first order transition between an interact- 
ing Fermi gas and a Fermi liquid, with its broad coexistence region extending from low 
densities up to about p - Signatures of this liquid-gas transition are nonetheless visible 
also in the chiral condensate at low temperatures and will be discussed in this work. 

Our starting point is the free energy density, F{p, T) = pF(p, T), of isospin-symmetric 
spin-saturated nuclear matter, with F(p, T) the free energy per particle. In the approach 
to nuclear matter based on in-medium chiral perturbation theory [7-9] the free energy 
density is given by a sum of convolution integrals of the form, 



oo oo 



P F (P, T )=4 J dppJC 1 n{p) + J dpi j dp 2 JC 2 n(p 1 )n(p 2 ) 
o oo 

oo oo oo 

+ J 'dpi J dp 2 J dp 3 IC 3 n(pi)n(p 2 )n(p 3 ) + pA(p,T) , (1) 



where K\,K, 2 and K 3 are one-body, two-body and three-body kernels, respectively. The 
last term, the so-called anomalous contribution A(p,T) is a special feature at finite 
temperatures [10] with no counterpart in the calculation of the ground state energy 
density at T = 0. As shown in ref.[7] the anomalous contribution arising in the present 
context from second-order pion exchange has actually very little influence on the equation 
of state of nuclear matter at moderate temperatures T < 50 MeV. 
The quantity 



p 2 /2M N ~ ■■- ' 
1 + cxp — 



(2) 
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denotes the density of nucleon states in momentum space. It is the product of the tem- 
perature dependent Fermi-Dirac distribution and a kinematical prefactor p/2ir 2 which 
has been included in n(p) for convenience. Mjv stands for the (free) nucleon mass. The 
particle density p is calculated as 

oo 

p = 4 J dppn(p) . (3) 
o 

This relation determines the dependence of the effective one-body chemical potential 
p(p,T; Mm) on the thermodynamical variables (p,T) and indirectly also on the nu- 
cleon mass Mat. The one-body kernel K\ in eq.(l) provides the contribution of the 
non- interacting nucleon gas to the free energy density and it reads [7]: 

K®= M " + *-£;-m' (4) 

The first term in K\ gives the leading contribution (density times nucleon rest mass) 
to the free energy density. The remaining terms account for (relativistically improved) 
kinetic energy corrections. 

The two- and three-body kernels, Ki and /C3, specifying all one- and two-pion exchange 
processes up to three loop order for the free energy density, have already been given in 
explicit form in refs.[7-9] and will not be repeated here. We recall from our earlier works 
that after fixing only a few contact terms the free energy density computed from these 
interaction kernels provides a realistic nuclear equation of state up to densities p < 2 pq. 
What matters in the following will be the dependence of the kernels Ki and /C3 on the 
light quark mass, m q , or equivalcntly, on the pion mass, m w , that is introduced by pion 
propagators and by pion loops. In-medium chiral effective field theory is the appropriate 
framework to quantity this pion-mass dependence in a systematic and reliable way. 

Application of the Feynman-Hcllmann theorem establishes an exact connection be- 
tween the temperature and density dependent in-medium quark condensate (qq)(p,T) 
and the derivative of the free energy density of (isospin-symmctric, spin-saturated) nu- 
clear matter with respect to the light quark mass m q . Using the Gell-Mann-Oakes-Renner 
relation m^f% = — m q (0\qq\0), one finds for the ratio of the in-medium to vacuum quark 
condensate 

(qq)(p,T) = P dF(p,T) 
(0\qq\0) ft dml ' { > 

where the derivative with respect to m 2 is to be taken at fixed p and T. The quantities 
(0|gg|0) (vacuum quark condensate) and f„ (pion decay constant) are to be understood 
as taken in the chiral limit, m q — > 0. Likewise, m 2 stands for the leading linear term in 
the quark mass expansion of the squared pion mass. 

In the one-body kernel Ki the quark (or pion) mass dependence is implicit via its 
dependence on the nucleon mass Mm- The condition dp/dM^ — applied to eq.(3) 
leads to the following dependence of the effective one-body chemical potential p on the 
nucleon mass Mm'- 



3 



^ - ' ip ^ = -4M N [ dp^- . (6) 



dM N 2M N Q% u J p 

o 

The nucleon sigma term <7n — (N\m q qq\N) = dMpf /dm 2 measures the variation of 
the nucleon mass with the quark (or pion) mass. Combining both relationships leads to 
the following m^-derivative of the one-body kernel: 

WCi = <tn_ ( 3p p 2 3p 4 \ 

dml mU + 2M N n% 3M 2 N 8M* J ' [ ' 

In the limit of zero temperature, T = 0, the terms in eq.(7) reproduce the linear decrease 
of the chiral condensate with density. The kinetic energy corrections account for the 
(small) difference between the scalar and the vector (i.e. baryon number) density. In the 
actual calculation we use the chiral expansion of the nucleon sigma term to order 
O(m^) as given in eq.(19) of rcf.[6]. The empirical value of the nucleon sigma term 
(at the physical pion mass = 135 MeV) is <jn = (45 ± 8) MeV [11]. Recent results 
for the quark mass dependence of baryon masses from lattice QCD and accurate chiral 
extrapolations [12] tend to give smaller values of this sigma term, but still consistent with 
the empirical <jn within errors. We use the central value <jn = 45 MeV in the present 
calculations, keeping in mind however that the existing error band in the determination 
of this quantity is still a primary source of uncertainties in following discussion. 

For the two- and three-body kernels, JC2 and /C3, related to one-pion exchange and 
iterated one-pion exchange, explicit expressions have been given in ref.[7]. Their deriva- 
tives with respect to the squared pion mass, dK,2^/dm 2 , are hence obvious and do not 
need to be written out here. The same applies to the anomalous contribution A(p, T) 
(see eqs.(14,15) in ref.[7]) and to the two- and three-body kernels related to 27r-exchange 
with excitation of virtual A(1232)-isobars (see section 6 in ref.[8]). In case of the one- 
pion exchange contribution we include the m T -dependent vertex correction factor T(m v ) 
as discussed in section 2.1 of ref.[6]. The short-distance contact term that produces a 
T-independent correction of order p 2 to the in-medium condensate is treated exactly in 
the same way as in ref.[6], i.e. all terms with a non-analytical quark-mass dependence 
generated by pion-loops are taken into account. 

Let us now turn to some three-loop contributions that are new and of special rele- 
vance for the in-medium chiral condensate. The first one comes from the pion self-energy 
diagram shown in Fig.l. It gives rise to the following m^-derivative of the two-body 
kernel: 



_ 3g\m\ 



dml ^UY 



{i 3 [(*r 2 ) 2 -(*i + 2 ) 2 ] 



+ (4£ 3 - 1) (X+ Xr 2 ) + (2*3 - 1) In || } , (8) 

with the abbreviations Xf- = [1 + (pi ± pj) 2 /m 2 ] -1 and the tttt low-energy constant 
£3 3. 

The chiral rnrNN contact vertex proportional to c\m\ generates 27r-exchange Hartree 
and Fock diagrams, also shown in Fig.l. Concerning the free energy density pF(p,T) or 
the equation of state of nuclear matter their contributions are actually almost negligible. 
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Fig. 1. Three- loop contributions to the free energy density of nuclear matter that are relevant for the 
in-medium chiral condensate. Upper diagram: pion self-energy correction; lower diagrams: two-pion ex- 
change Hartree and Fock terms involving the ttitNN contact interaction proportional to the low-energy 
constant c\. 




Fig. 2. Density dependence of the chiral condensate in isospin-symmetric nuclear matter at temperature 
T = 100 MeV. Starting from the linear density dependence (dashed curve) characteristic of the free 
nucleon Fermi gas, the following interaction contributions are successively added: one-pion exchange 
Fock term (lw), second order (iterated) pion exchange, irreducible two-pion exchange (2-7r), two- and 
three-body contributions from 27r exchange with intermediate A excitations (A), and two-pion exchange 
with ttttNN vertex involving the low-energy constant ci. Pauli blocking effects are included throughout. 

However, when taking the derivative with respect to m^. as required for the calculation 
of the in-mcdium condensate, these contributions turn out to be of similar importance as 
other interaction terms. The corresponding contribution to the derivative of the two-body 
kernel reads: 



dml Stt/* 
with the auxiliary function 
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G(x) = 8x(3 + x 2 ) arctanx - 5 ln(l + x 2 ) - 100x 2 . (10) 

The 27r-exchange Hartree diagram with one cimj-vertex contributes the following piece 
to the three-body kernel: 

dKtfL = _ x _ )(Xi+2 + x - - 3) + In f|} , (11) 

while the three-body term associated with the 27r-exchange Fock diagram with one c\m 2 - 
vertex gives: 



pi + g§ ^_jm| ln p + (pi + p3)x + + (pi _ p3)x _ 



p 3 4p| X+_ 

„2 ^2 q^,2 

(12) 

L 13 J 

Last not least we incorporate the effects of thermal pions. Through its m^-derivative the 
pressure (or free energy density) of thermal pions gives rise to a further reduction of the 
T-dependent in-medium condensate. In the two-loop approximation of chiral perturba- 
tion theory including effects from the 7T7r-interaction one finds the following shift of the 
condensate ratio in the presence of the pionic heat bath [13-15]: 



S(qq)(T) 3ml 



(0\qq\0) (2ttA) 2 



(13) 



with the functions Hi^m^/T) defined by integrals over the Bose distribution of thermal 
pions: 

oo oo 

HUy) = [ dx 1 , H 3 (y) = y~ 2 [ dx ^ ~f . (14) 

v ' y 

We proceed with a presentation of results. As input we consistently use the same 
parameters in the chiral limit as in our previous works [6], namely: f n = 86.5 MeV, 
g A = 1.224, ci = -0.93 GeV^ 1 and M N = 882 MeV. Concerning the contact term 
representing unresolved short-distance dynamics, we recall from ref. [6] that its quark 
mass dependence, estimated from recent lattice QCD results [16], is negligibly small 
compared to that of the intermediate and long range (pion-exchange driven) pieces. 

It is worth pointing out again that in-medium chiral perturbation theory with this 
input produces a realistic nuclear equation of state [8,9], including a proper description 
of the thermodynamics of the liquid-gas phase transition. Apart from temperature T, 
the additional "small" parameter in this approach is the nucleon Fermi momentum pp in 
comparison with the chiral scale, 4^/^ ~ 1 GeV. Our three-loop calculation of the free 
energy density is reliable up to about twice the density of normal nuclear matter. It can 
be trusted over a temperature range (up to T ~ 100 MeV) in which the hot and dense 
matter still remains well inside the phase of spontaneously broken chiral symmetry. 

Fig. 2 shows a representative example, at T = 100 MeV, displaying stepwise the effects 
of interaction contributions to the density dependence of (qq) (p, T) arising from the 
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m^-derivative of the chiral two- and three-body kernels K.2,3- As in the T = case 
studied previously [6], the pion-mass dependence of correlations involving virtual A (1232) 
excitations turns out to be specifically important in delaying the tendency towards chiral 
symmetry restoration as the density increases. Once all one- and two-pion exchange 
processes contributing to dK,2/dm^ and dJC^/dm^ are added up, the chiral condensate 
at T = 100 MeV recovers the linear density dependence characteristic of a free Fermi gas. 
However, this recovery is the result of a subtle balance between attractive and repulsive 
correlations and their detailed pion-mass dependences. Had we taken into account only 
iterated one-pion and irreducible two-pion exchanges, the system would have become 
instable not far above normal nuclear matter density as can be seen in Fig. 2. In fact, this 
instability would have appeared at even much lower densities in the chiral limit (m„ — ¥ 
0). This emphasizes once more not only the importance of terms including A(1232)- 
excitations, but also the significance of explicit chiral symmetry breaking by small but 
non-zero quark masses in QCD and the resulting physical pion mass, m v = 135 MeV, in 
governing nuclear scales. 




Fig. 3. Ratio of chiral condensate relative to its vacuum value as function of baryon densityp in 
isospin-symmetric nuclear matter, for different temperatures up to T = 100 MeV as indicated. The 
effects of thermal pions are included. The dashed line at T ~ results through a Maxwell construction 
in the coexistence region of the nuclear liquid and gas phases. 

Fig. 3 shows the systematics in the variation of the chiral condensate with temperature 
T and baryon density p. These results include all nuclear correlation effects and also the 
(small) additional shift from thermal pions. The latter correction is visible only at the 
highest temperature considered here (T = 100 MeV) where the chiral condensate at zero 
density starts to deviate from its vacuum value. The actual crossover transition at which 
the condensate drops continuously to zero is around T ~ 170 MeV [17]. . 

At zero temperature, the hindrance of the dropping condensate at densities beyond 
normal nuclear matter comes primarily from three-body correlations through K3 which 
grow rapidly and faster than IC2 as the density increases. The heating of the system 
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reduces the influence of £3 relative to K.2 continuously as the temperature rises, so that 
their balance at T = 100 McV produces a small net effect in comparison with the free 
Fermi gas. 

At T = 0, the solid line in Fig. 3 does not yet take into account the fact that the density 
range up to and including normal nuclear matter density covers the coexistence region of 
the nuclear liquid and gas phases [9] . Any first-order phase transition is expected to leave 
its mark also in other order parameters, and so it does for the chiral quark condensate. 
Based on the usual Maxwell construction, the dashed line in Fig. 3 indicates this effect. 
It becomes much more pronounced when the chiral condensate at low temperatures is 
plotted as a function of the baryon chemical potential: 

H = M N +(l + p^F(p,T). (15) 

The discontinuity indicating the first-order liquid-gas transition at T smaller than the 
critical temperature for this transition, T c ~ 15 MeV, is clearly visible in Fig. 4. At this 
point our results are consistent with a recent investigation aimed at an understanding of 
chemical freeze-out in heavy-ion collisions at large baryon densities [18], where a similar 
effect is found using a chiral meson-baryon model Lagrangian. Another side effect induced 
by the first-order liquid-gas phase transition in the discussion of the in-medium chiral 
condensate (qq)(p) is that the frequently advocated "low-density theorem" needs to be 
modified: 

(0\qq\0) miff' ( b) 

where the effective nucleon sigma term ctjv — 36 MeV measures the quark mass depen- 
dence of the sum Mm + Eq, with Eo ~ —16 MeV the binding energy per particle of 
saturated nuclear matter. The usual version of eq.(16) with the nucleon sigma term 
in vacuum assumes that at sufficiently low densities nuclear matter could be treated as 
a non-interacting Fermi gas. 

In summary, this is the first calculation of the quark condensate at finite temperature 
and density that systematically incorporates chiral two-pion exchange interactions in the 
nuclear medium. Correlations involving intermediate A(1232)-isobar excitations (i.e. the 
strong spin-isospin polarizability of the nucleon) together with Pauli- blocking effects 
are demonstrated to play a crucial role in stabilizing the condensate at densities beyond 
that of the nuclear matter ground state. The results reported here set important nuclear 
physics constraints for the QCD equation of state at baryon densities and temperatures 
that are of interest e.g. in relativistic heavy-ion collisions. In particular, we find no indi- 
cation of a first-order chiral phase transition at temperatures T < 100 MeV and baryon 
densities at least up to about twice the density of normal nuclear matter. 
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Fig. 4. Ratio of chiral condensate relative to its vacuum value as function of baryon chemical potential in 
symmetric nuclear matter at low temperatures characteristic of the nuclear liquid-gas phase coexistence 
region. 
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